On 
ON 



(N 



> 



CU 



43 



X 



July, 1994 



Modified tetrahedron equations and related 3D integrable models 



H.E. Boos,V.V. Mangazeev,S.M. Sergeev 

Institute for High Energy Physics, 142284 Protvino, Moscow Region, Russia 



o, 

Abstract 



Using a modified version of the tetrahedron equations we construct a new 
family of iV-state three-dimensional integrable models with commuting two- 



layer transfer-matrices. We investigate a particular class of solutions to these 



equations and parameterize them in terms of elliptic functions. The corre- 
sponding models contain one free parameter k - an elliptic modulus. 



1. Introduction 

Recent progress in the construction of solvable lattice models in higher di- 
mensions has stimulated efforts to find new interesting examples "living" 
in D > 2 dimensions. As is known the main ingredient of practically all 
two-dimensional solvable models is a Yang-Baxter equation (YBE) which 
provides an existence of a family of commuting transfer-matrices. As a rule, 
the presence of such a family permits us to calculate the partition function 
per site in the thermodynamical limit. Unlike two-dimensional conformal 
field theories the theory of the YBE can be easily generalized for the cases 
D > 2 [1-5] and it will lead us to the tetrahedron equations (D = 3) and 
.D-simplex equations. The main problem, of course, is to solve them, because 
the number of unknowns and equations grows with increasing dimension in 
a catastrophic way. 

Recently Bazhanov and Baxter in their paper [7] have generalized the 
trigonometric Zamolodchikov model [1-3] for the case of N colours. The 
weight functions of their model satisfy tetrahedron equations ([8]) and pos- 
sess some remarkable symmetry properties. But the essential shortcoming of 
Bazhanov-Baxter iV-color model is an absence any temperature-like param- 
eter. 

Authors of [9, 14] have proposed a new class of three-dimensional models 
containing a temperature-like parameter k. In the limit k — > the weight 
functions of these models reduce to the "static limit" case of the Bazhanov- 
Baxter model. The main idea of papers [9, 14] is to introduce a couple of mod- 
ified tetrahedron equations (MTE). In two dimensions this idea have been 
considered in [6] . Then using a couple of modified tetrahedron equations, the 
authors of [9, 14] have constructed two-layer commuting transfer- matrices for 
the checkerboard lattice of weights. 

This paper, in fact, is a continuation of [9,14]. Our main purpose is to 
remove the "static limit" condition which is a constraint on three "angle" 
parameters for weight functions. The result is that each elementary weight 
function of the model will depend on three spectral parameters and one 
"modulus". Unfortunately, to remove the "static limit" condition we were 
forced to refuse from the second modified tetrahedron equation. It leads to 
a more complicated commuting family of two-lawer transfer-matrices. 

The paper is organized as follows: In Section 2 we give necessary defini- 
tions, the explicit expression for weight functions and write out the modified 



tetrahedron equations. In Section 3 we recall the proof of the modified tetra- 
hedron equations based on the Star-Square relation, give all necessary alge- 
braic constraints on the parameters and consider an important submanifold 
of "equal" moduli in the space of parameters (in next sections we restrict 
ourselves to this case only). In Section 4 we introduce a parameterization of 
weight functions in terms of angle-like variables. In fact, all the formulae of 
Section 4 can be rewritten in terms of elliptic functions and we give the cor- 
responding formulae in the Appendix. Section 5 contains a discussion of the 
symmetry properties of the weight functions under the group of transforma- 
tions of the elementary cube. In Section 6 we construct a commuting family 
of two-layer transfer-matrices using the solutions obtained from the modified 
tetrahedron equations and discuss the structure of the lattice model. Finally, 
the Appendix contains all the elliptic formulae and a detailed consideration 
of the N = 2 case. 



2. Modified tetrahedron equations and Body- 
Centered-Cube (BCC) ansatz for weight 
functions 

In this section we present the modified tetrahedron equations (MTE) and 
briefly discuss their properties. We recall necessary definitions from [8, 11, 14] 
and give the symmetrical form of the weight functions W 

Further we will follow a statistical mechanics interpretation and use the 
notations introduced by Baxter in [3] for the Zamolodchikov model. Namely, 
consider a simple cubic lattice C and at each site of C place a spin variable 
taking its values in Z^, for any integer N > 2. 

To each elementary cube of the lattice we can assign some weight function 
W(a\efg\bcd\h), depending on eight surrounding spin variables (see Fig. 1). 




W(a\e,f,g\b,c,d\h) 



Fig. 1 

Note that for different elementary cubes of the lattice L we can use differ- 
ent weight functions W^\ W^ etc. Later we will fix the explicit structure 
of the lattice C, demanding that the two-layer transfer-matrices commute 
between themselves. 

Now let us consider two sets of weight functions: W, W, W", W" and 
W, W , W , W . Suppose that they satisfy the following equations: 

Y^ W(a 4 \c 2 , c u c 3 \bi, b 3 , b 2 \d)W (ci|fe 2 , a 3 , 6i|c 4 , d, c 6 |& 4 ) 

d 

x W"(bi\d, c 4 ,c 3 \a 2 ,b 3 ,b 4 \c 5 )W (d\b 2 , h, b 3 \c 5 , c 2 , c 6 |oi) 
= H W'"(bi\ci, c 4 , c 3 |a 2 , a 4 , a 3 |<i)W^ (ci|fe 2 , o 3 , a 4 |rf, c 2 , c 6 |oi) 



x 



W^'(o4|c 2 , d, c 3 |a 2 , 6 3 , ai|c 5 )W(d|ai, a 3 , a 2 |c 4 , c 5 , c 6 |6 4 ), (2.1) 



where a iy bi, Cj G Zjy, i = 1, . . . , 4, j = 1, . . . , 6. We will call relations (2.1) 
the modified tetrahedron equations. Note that all eight weights in (2.1) 
are independent, but we placed W and W in relations (2.1) so as to follow 
the notation of papers [9,14]. For the case W = W we come to the usual 
tetrahedorn equations (see, for example, [3]). We would like to stress that 
contrary to [9, 14] we do not demand the validity of the dual variant of (2.1) 
(with all IV's replaced by W's and vice versa). The absence of the dual 
variant of (2.1) permits us to overcome the static limit condition (see [9, 14]) 
but leads to a more complicated commutativity of the two-layer transfer- 
matrices. 

Let us recall some definitions used in papers [8, 11, 14]. First denote 



u = exp(2ni/N), u 1/2 = exp(ni/N). 



(2.2) 



Further, taking x, y, z to be complex parameters constrained by the Fermat 



equation 

and / to be an element of Z^, define 

w (x,y,z\l) = l[—^. (2.4) 

In addition, define the function with one more argument 

w(x, y, z\k, I) — w(x, y, z\k — l)&(l), k,l E Z N , (2.5) 

where 

$(/) = J^+ N )l\ (2.6) 

Let us mention also two formulae for the w functions, which are useful for 
calculations: 

w(x, y, z\l + k) = w(x, y, z\k)w(xu k , y, z\l), (2.7) 

w(x, y, z\k, I) = u kl /w(z, u 1/2 y, ux\l, k). (2.8) 

Now introduce the set of homogeneous variables Xi, % — 1, . . . , 8 and xi 3 , 

X2i, x 58 , x &7 satisfying 

T N _ N _ N N _ N _ N N _ N _ N N _ N _ N (l) q n 

We will also need six additional variables to ensure the necessary transforma- 
tion properties of the weight functions under the group G of transformations 
of a three-dimensional cube. Define them as 

U N = (x 3 X 5 /xi) N — Xg , V N = Xj — (X4X 6 /X2) N , 

£ N = ( Xl x 7 /x 3 ) N - x%, X N = xf - (x 2 x 8 /x 4 ) N , 

fJ> N = (x 13 X 24 X 7 /x 3 X 4 ) N - (x 58 X 67 /x 8 ) N , 

V N = (Xl3X24X 6 /x 1 X 2 ) N - (x 58 X 67 /x 5 ) N . (2.10) 

Using these definitions we define the weight function W(a\efg\bcd\h) as 

~w(x 58 x 67 , x 8 fx, x 13 x 24 x 7 x 8 /x 3 x 4 \a + d,e + /)l 1 I 2 



W(a\eJ,g\b,c,d\h) 



w(x 58 x 67 , x 5 v, x 13 X24X 5 x 6 /x 1 x 2 \g + h,b + c) 



X 



w(xix 8 , Xiu, x 3 x 5 \e + h,d+ c) 



w(x4X 6 , x 2 v, x 2 x 7 \a + b,f + g) 



1/2 



w (x 2 x s , x 4 A, x^x 5 \e + g,a + c) 



w(a; 3 a; 6 , a; 3 £, Zi^jfr + d, f + h) 



1/2 



u bf u (ag+gin-bh)/2 ^ w(x 3 , x 13 , x^d, h + a)w(x 4 , x 24 , x 2 \a, g + a) 



u ag u (dh+de+ea)/2 l^itf^g, X 58 , X 5 |e, C + a)w(x 7 /cc>, X 67 , X 6 |/, 6 + (j) J 

(2.11) 

where the lower index "0" after the curly brackets implies that the expres- 
sion in the curly brackets is divided by itself with all exterior spin variables 
equated to zero. 

Note that all gauge multipliers (face and edge types) before curly brackets 
in (2.11) are necessary in order to provide the correct symmetry properties 
of the weight functions under the transformations of the group G. 

Formula (2.11) generalizes the weight functions of the model proposed by 
Bazhanov and Baxter [7, 11]. In fact, it coincides with the weight functions 
from [14] up to gauge multipliers and we will also call (2.11) the Body- 
Centered-Cube (BCC) ansatz for weight functions. But in contrast to [14] 
we will introduce a new parameterization for xi in a such way that each 
weight function will depend on three independent spectral parameters and 
modulus parameter k. 

3. The proof of the modified tetrahedron equa- 
tions 

In this section we recall the proof of the modified tetrahedron equations (2.1) 
for weight functions (2.11). In fact, it reproduces the method of papers [8, 14] 
and we give it here for completeness. 

The main idea of papers [8, 14] is to reduce complicated modified tetra- 
hedron equations (2.1) to a couple of much more simple relations: namely 
"inversion" and Star-Square ones. 

The "inversion" relation for functions w(x, y, z\k, I) has the form: 

^ w(x,y,z\k,l) _ (1 ~ z/x) 

where I, m E Zn, x,y,z satisfy (2.3) and 5^ m is the Kronecker symbol on 

Zn- 
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The Star-Square relation permits us to calculate the sum over the one spin 
variable from the product of four w functions provided that their arguments 
satisfy some algebraic constraint. To avoid the introduction of additional 
variables for ensuring the cyclic property of all u>'s modulo N, we will use a 
non-cyclic analog of w function, defined recurrently as follows: 

"(*|0) = i, Tlf '°n = 7T J — IT' /GZ ' (3 ' 2) 

w[x\l — 1) (1 — XUJ 1 ) 



where Z is the set of all integers. It is obvious that 

w(x, y,z\l) — (-) w(x/z\l), I G Z N , (3.3) 



where index /, being considered modulo N, is interpreted as an element of 
Z N . 

Then the Star-Square relation can be written as: 

y, w(x u g/i, zi\a + a)w(x 2 , y 2 , z 2 \b + a) } 

<£z N w ( x 3' Vzi z z\ c + cr)w(x 4 , 2/4, z A \d + a) Jo 

(x 2 z 1 /x 1 Z2) a ~ b (x 1 y 2 /x2Z 1 ) b (z 3 /y 3 ) c (z4/y 4 ) d 

$(a - b)uj l ~ a+b y 2 

w(ux 3 x i ziz 2 /xix 2 z 3 z i \c + d — a — b) a a\ 

x —TiZYT " x 3 Z 2 - a-Vi ^2,, 7T ' ^ A ) 
w( 1 a — a)w{ \c — b)w{ \c — a)w{ \d — b) 

X\Z^ X 2 Z 3 X\Z 3 X 2 Z^ 

where the lower index "0" after the curly brackets indicates that the l.h.s. of 
(3.4) is normalized to unity at zero exterior spins, and the following constraint 

is imposed 

y\ 2/2 z 3 z 4 . s 
= uj. (3.5) 

Z\ Z 2 J/3 J/4 

Note that the separate u>'s in the r.h.s. of (3.4) are not single- valued functions 
on Zn, while the whole expression is cyclic in the exterior spins a, b, c, d. 

Relations (3.1), (3.4) are proved by using properties of Fourier transfor- 
mation over spin variables and the detailed proof was given in [8] . 

Now we will prove relations (2.1) for ansatz (2.11) with some appropriate 
choice of constant multiplier R. Instead of weights W, W, W", W" let 
us substitute into (2.1) explicit formula (2.11) with corresponding sets of 



parameters: 

w, w, w", w" -> w( Xi , x tJ ),w(x' t , 4), w(4, 4), ^(4, 4), 

W 7 , W,W',W" - ^(x,,%), w(^,4), w(^',4), w(^",4). (3.6) 



Now we will show that all equations (2.1) are equivalent to some algebraic 
system of nonlinear equations for parameters Xi, Xi etc. 

Hereafter we will fix a normalization of all parameters x's that enter in 
the definition (2.11) as 

x 3 = 1, x 4 = 1, :r 7 = 1, x 8 = 1, (3.7) 

and similarly for all the x's. 

It is easy to see that after substituting of (2.11) in relations (2.1) we will 
explicitly obtain 24 pairs of w functions coming from the multipliers before 
the curly brackets in (2.11). Let us demand that all these factors cancel each 
other (for each w function there is one partner with the same dependence 
from spin variables). Then we obtain 24 relations for variables Xi, Xij and 
x's, x'ijS (only 21 of the constraints are independent) and 24 relations for the 
additional variables defined by (2.10). In fact, they relate the set of rapidity 
like parameters between different weight functions. 

Let us multiply both sides of (2.1) by the following product of w weights: 



w (x" , a/ fi ' 7 , a/ fi ' | c 4 , a 2 + h ) w(x' 7 , ufx^ , x' 6 



w(x , z/uj, x'I^x'I'Icq, ai + /i) w(x4, x 24 , ux 2 



03,C 4 + /3 



a 4 , C3 + Z 4 ) 



x ^(xs xga xlju)\b 2 , c 2 + m 2 ) w{ux' H , x^, a^|c 2 , 6 3 + m 3 ) ,3 g x 

w(xl', X24, oux'^'lh, c 3 + mi) w(x 3 /u, x 13 , £i|c 6 , 64 + m 4 ) l ' ; 

and sum over 01,02,03, a 4 , 01,62,63,64- Note that due to the "inversion" 
relation (3.1) we do not lose any information after such transition of spins 

Oj, h to /j, TOj. 

The functions w in expression (3.7) are chosen in such a way that using 
relation (3.1) we can calculate the sums over the spins 01,02,03,04 in the 
l.h.s. and those over 61,62,63,64 in the r.h.s. of the obtained equation and 
cancel the summations over the spins ex's, which come from expression (2.11) 
for the functions W's and W's. Now let us consider the applicability con- 
ditions (3.5) of Star-Square relation (3.4) for the sums over ai,a 2 , 03, a 4 in 
the r.h.s. and over 61,62,63,64 in the l.h.s. of the obtained equation. We 



obtain eight conditions on the x's and x's. Applying relation (3.4) eight 
times and calculating the sums over d spin in the l.h.s. and r.h.s (the spin 
structure of the sums over d has the form of relation (3.1) and we demand 
that corresponding variables x, y, z entering in the arguments of functions 
w(x, y, z\k, I) are constrained in such a way that relation (3.1) can be used), 
we come to the equation without any summation. The l.h.s. and r.h.s. of 
this equation consists of the products of ten w functions with the same spin 
structure and expressions like x™ 2 coming from relation (3.4). Let us impose 
all the necessary constraints on parameters Xi, Xij, Xi and x~ij to satisfy this 
equation. The explicit calculations show that we obtain 11 more indepen- 
dent constraints on parameters. Taking into account 21+24=45 constraints 
coming from a cancellation of "face" factors in (2.11) we have an algebraic 
system of 56 nonlinear equations. 

Taking into account relations (2.9,2.10) it is easy to understand that 
only 23 of the constraints are independent on the level of iV-th powers of 
coordinates x's and x's. Since each weight of form (2.11) depends only on 
four independent parameters, we have a nine parameter solution of relations 

(2-1). 

To order all these constraints let us introduce the following rapidity like 
parameters: 

i 3 = a , (3.9) 

^58^67 

^i. (3.10) 

^58^67 

and the same for the x's. Then the set of 24 relations following from a can- 
cellation of "face" factors (which do not contain additional variables defined 
in (2.10)) has the form: 

(3.11) 

(3.12) 



(3.13) 



X 5 X 6 

xix 2 : 


k = 


Xq 
X 2 


12 = 


_ x 5 

X2 


h = 


Xi 

■-UJ-, 

x 5 


h = 




h : 



h = %', 


*2 = i'i, 


h = i'2, 


l l — l 3 i 


L 3 ~ L 3i 


A = 


ji'i 
n i 


— —in 
3i = 32 i 


32 = 32i 


— —a 

33 = 32i 


—I —in 

3\ = 33 ' 


—l —ll 

33 = 33i 


—a 
Ji 


—in 
= Ji- 


-in 

3l=*2i 


h = H, 


-II 
33 = hi 


3i = hi 


., -II 

33 — Hi 


i'i 


-in 


h = j'2, 


12 = f 2 , 


~ -II 
*3 = 32 > 


-i -in 
h=3 3 > 


-1 -II 

h — 33i 


-n 

h 


- i'" 

— 3\ 



Only 21 of the 24 relations (3.11-3.13) are independent. A simple consequence 
of (3.11-3.12) is 

a a a" a'" ,.,., 

A = i = ? = i* = ^ (3 - 14) 

and we can conclude that a parameter A is an absolute invariant of relations 
(2.1). 

We will also define rapidity variables of the mixed type: 

^13^67 ^13^58 X\ , 

«i = — ^, a 2 = — ^, «3 = — ■ (3.15) 

^67^13 ^58^13 x l 

Then six of the eleven additional relations coming from the coincidence of 
the structures of l.h.s. and r.h.s. of the modified tetrahedron equations (after 
using of "inversion" and Star-Square relations) take the form 

/// / // / /// / // // III /o 1 r\ 

a,\ = a 2 , a 2 = Oi 2 , 0:3 = a 2 , a 1 = a 3 , a 3 = a 3 , a 1 — a 1 . (o. lb) 

Now we can similarly rewrite 24 of the constraints containing additional 
variables defined in (2.10). Define 



01 = 


1/2 X n UX 5 
U) * V, f3 2 = , 03 = 

X2 ^58^67 




(3.17) 




1/2 ux 1 £ \x 




(3.18) 


7i — 


<*> ? 12 — j 73 — 

X5 X6 ^58^67 


and similarly for f3 { , 7 y i . 


Then we have 






0i = 2 , 02 = 2 , 


03 = 0^', #=03", 03 = 03, 01' = 


= 0i", 


(3.19) 


— — in — —i 

7i = 7 2 , 7 2 = 7 2 , 


_ _// — 1 —in — 1 — 11 — n 

7 3 = 7 2 , 7i = 7 3 , 7 3 = 7 3 , 7i = 


= 7? 


(3.20) 


and 








15'" 73' 

7i = P 2 , 72 = 2 , 


73 = 0'2, 7l=0 3 ", 73 = 03, H- 


= 0i', 




01=72', 02=72, 


03 = 7^', 0;=73", % = ll 01 


= 7i"- 


(3.21) 



The remaining five constraints include parameters of all four sets of weight 
functions and correspond to a generalization of tetrahedron quadrilateral 
constraint for the Zamolodchikov model. We can choose them as 

X 2A X 2 4 : = X 2 4 X 2 A \O.ZZj 
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and 

~ H7 — // — /// ' ^^7 77 777" ' ^ — 77~~ ' 7 — w ■*■■ 

Xi X 13 J/^ X 13 X 2 4 X 2 X 24 X24 -^67 x 58 J '24 

(3.23) 
Now let us briefly discuss the main properties of relations (3.11-3.23). 
As we mention all these equations have a nine parametric solution. We can 
construct it in the following way. Relations (3.11,3.22) give 7 restrictions on 
16 parameters x iy x\, x" and x'" . Then all parameters xi can be calculated in 
a unique way using relations (3.12,3.13,3.16,3.22). All remaining constraints 
will be satisfied automatically up to some JV-th root of unity. And choosing 
the degrees of the cu's in a corresponding way we will satisfy all relations 
(3.11-3.23) and therefore the modified tetrahedron equation (2.1). 

Now let us consider an important particular case of relations (3.11-3.23): 

a = a ' = a" = a". (3.24) 

Constraints (3.24) reduce the number of parameters to six. In this paper we 
will consider only a class of solutions for the modified tetrahedron equations 
restricted by conditions (3.24). Then the parameter a can be interpreted 
as an invariant. Detailed analysis of equations (3.11,3.22) shows that on 
the surface (3.24) we have two solutions. Only the first one contains the 
Zamolodchikov - Bazhanov - Baxter model in the limit a — > 1 . So we will 
investigate this case only. Then one can show that 

A = a 2 = =2. (3.25) 

a 

It is easy also to show that the following relations between the iV-th powers 
of the parameters ik, jk, otk, Pk and 7^ are valid: 



i* j N ■ 

3k = Tap a k = ~n 7' Pk = h ~ 1j Ik = Jk ~ 1) k — 1,2,3. 



7 iV' k „-iV _ 

(3.26) 

Note that formulae (3.26) for the iV-th powers of parameters a^ give a con- 
venient way to find all parameters x's in terms of x's. 

In the next section we will give a convenient parameterization of all pa- 
rameters in terms of variables resembling tetrahedron angles. 
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4. Parameterization 

Define m, T\, T 2 , T 3 through the relations 

-mT? 



x 2 
X* 



%2 



-mTl 



±13±24 
x 58^67 



T N N 
_ mT 2. £l^2. 



3 ! 



X 5 X 6 



m 



(4-1) 



The advantage of this parameterization is that when we put m = 1 and Tj = 
tan(0j/2), we obtain the conventional parameterization of the Zamolodchikov 
- Bazhanov - Baxter model. 

Solving (4.1) with respect to the xf 's, we obtain 



x 



N 



1 



mTl 



T i T 2 \ 1 + m _1 7T: 



1-2 exp(m 3 J 



.7; 



;r 



iV 



iV 



T\T 2 

-iT, 



1 + mTi 



\l + wT x T, 



2 exp^a 3 ; 






iV 



-m 



-m 



"•#3 



T 2 \ 1 + m -1 ^ 



2 exp(m 3 J 



-iT 2 



1 + mT' 



t i ^ 1 + m _1 r 3 

where a 3 is defined from its cosine: 



1-2 exp^a 3 j 



(4.2) 



cos(a 3 ) 



1 + T 2 T 2 2 - T 2 T 2 - T 2 2 T 2 



2T 1 T 2 J(l + mTi)(l+m- l Ti 



(4.3) 



The formulae for the x^s can be obtained from these ones by changing m — : 
to -1 , a 3 being the same. 

The formula (4.3) can be rewritten in another form by the substitution 



Ti 



Up ~ ^p 



K J|) 



Up + C p 



(4.4) 
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where 

D p — 1 — K D p , C p = 1 — D p , 

m= l^|, k 2 + k' 2 = l. (4.5) 

and p — 1, 2, 3. Then we have (defining the angles Oi and a2 at once) 



COS O r 



fc bpbq 



(4.6) 



where {p, g, r} = {1,2, 3}. Introduce else the angle a as follows: 

D 1 D 2 D 3 + k 2 C x C 2 C z 



cos(a ) 



k 



a 



(4.7) 



(4i 



ao is symmetric under the permutation of T\,T 2 ,T 3 , and the sign of ao we 
choose as follows 

sin(a ) = ksm(a r )S p S q . 

To write down x^ we need to introduce a "linear excesses" 

n oi + a 2 + a 3 + a a„ + a„ - a r - a 
A) = vr , /3 r = -^ , 

— oi + a 2 + a 3 — a -^ a p + a q — a r + a 

P0 = 71 o ' ^r = o • 



(4.9) 



Then 



x 



N 
13 



;r 



-m 



1/2 ?]. 



T^\ 



1t^2\ 



l + m" i T 1 2 )(l+m" 1 T ; 



;i+^ 1 t 3 2 ) 



■exp(-iA)), 



iV 

21 



-m 



-1/2T, 



l + mT 2 )(l + mT 2 2 



\ (1 + m^TT 



exp(-i/? 3 ) 



*£ = "^~ 1/2 ^ 



\ 



(l+m^Tl) 



exp(i/3 2 ), 



AT _ _ rn -l/2 1 

7Y\ 



Xg 7 = — m 



;i+mT/)(l+m- i T 2 / ) 
(l+m-^ 2 ) 



exp^), (4.10) 
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We mention again that the xf's and x^s can be obtained by the same 
formulae as xf and xfj with m replaced by mr x (hence k replaced by —k 
and a replaced by — a ), ^1,2, 3 being unchanged. 

One can take the iV-th roots of the expressions (4.2,4.10) so that when 
m — ► 1 there appears the parameterization of the Zamolodchikov - Bazhanov 
- Baxter model. Namely, if m, T r and a r , r = 1,2,3, are real positive then 
we can choose the phases as follows 

x\ = \xi\ exp(ia 3 /N), X13 = Ixizlcu 1 / 2 exp(— i/3 /N); 
x 2 = \x 2 \exp(ia 3 /N), x 24 = \x24\oj 1/2 exp(-i(3 3 /N); 
x 5 = \x 5 \u 1/2 exp(ia 3 /N), x 5S = \x 58 \uj 1/2 exp(i(3 2 /N); 
Xq = \x 6 \uj~' 1 / 2 exp(ia 3 /N), x 6 j = |x 67 |o; _1 / 2 exp(i/3 1 /A^). (4.11) 

In addition we reproduce some exotic formulae for the "spherical sides" . 



sin /3q sin fa sin /3 2 sin f3 3 

m, 



\ sin P sin fi x sin (3 2 sin (5 3 
mT 2 = sin/3psin/j g 
sin /3 sin (3 r ' 



= smapsma^ 
sin a p sin a g 

Now rewrite the scheme of arguments of weights (W, W) in the modified 
tetrahedron equations. Exhibiting the dependence on T 1]23 and correspond- 
ing a ,i, 2 ,3 as 

W = W(T U T 2 , T 3 ) -> a , 01, a 2 , a 3 , (4.13) 

the arguments of the weights in the modified tetrahedron equations can be 
written as follows 

W = W(T U T 2 , T 3 ) -»■ a , ai, a 2 , a 3 ; 
W" = W(T 4 , T 2 , T 5 ) -> o'o, a[, a' 2 , a 3 ; 
W" = W(T 6 , 1/T 3 , T 5 ) - <, a'/, a 2 ', a 3 '; 
W" = W(T 6 , T l5 1/T 4 ) - o' ", a'/', <, a 3 ". (4.14) 

In terms of the "spherical sides" of the "tetrahedron" , the modified tetrahe- 
dron equations correspond to the set of equivalent relations 

ao + a 2 + a 2 " = a'2, a' + a 2 + a 3 = o 3 , 
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a'o + 03 + a?" = o' 3 , a'o + en + a" = a[ (4.15) 

Concluding this section we note that some of the formulae resemble those for 
the elliptic functions. Surely one can put S = sn(#, k) for S and k defined 
by (4.4,4.5). Such version of the formulae we shall collect in the Appendix, 
where we also give the explicit form of the Boltzmann weight functions for 
the case N = 2. 



5. Symmetry properties 

In this section we give all the necessary formulae for the transformation 
of parameters from which the weight function (2.11) depends on under the 
action of elements of the group G of all symmetry transformations of the 
three-dimensional cube. The main advantage of introducing parameters (2.9- 
2.10) is that any element from group G induces a multiplicative law of the 
transformation of parameters x iy Xij, u, v, £, A, p and v. 

Any element of the group G can be represented as a composition of two 
generating elements: r and p (see [11, 13]). The action of these two elements 
on the set of spins {a\e, f, g\b, c, d\h} (see Fig. 1) can be described as follows 

r{a\e, f,g\b, c, d\h} = {a\f, e, g\c, b, d\h} (5.1) 

and 

p{a\e, f, g\b, c, d\h] = {g\c, a, b\f, h, e, \d}. (5.2) 

In fact, r corresponds to the reflection of a three-dimensional cube in the 
plane {aghd} and p - to a rotation of tt/2 around the vertical axis of the 
cube. Further it will be convenient to remove normalization conditions (3.10) 
and to restore the homogeneity of the parameters x's. 

The action of the element p can be obtained by use of the Fourier trans- 
formation (see [11]) and the action of r is a quite obvious. Here we give 
only the resulting formulae for their action on the coordinates. Note that 
all coordinates can be split into four multiplets: (x3,xi3,xi), (x4,X24,a;2), 
(x$,X5s,X5,u, A) and (x 7 ,x Q7 ,xq,v ,£,/i, v) and that any group element acts 
homogeneously on these four sets. Also, we must demand a trivial action of 
the elements: r 2 , p 4 and (rp) 6 (trivial transformations) on all coordinates. 
This will fix some uncertainties in the iV-th roots of unity in the transforma- 
tion laws. 
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Then a correct action of the elements r and p looks like: 

(x 3 , x 13 , Xi) -H- (x 3 , x 13 , xi), (x 4 , x 24 , x 2 ) A (x 4 , x 24 , £2), 
(x 8 , x 58 , x 5 , w, A) A (x 7 /w, a; 67 , s 6 , uj~ 1/2 ^, u~ 1/2 ^), (5.3) 

(aj 7 , x 67 , x 6 , v, £, fi, v) A {ux % , x 58 , x 5 , ^ 1/2 ^-, ^ 1/2 ^p w^, ^f ) 
and 

(X 3 , X13, Xi) A (x 8 Xi3, XiU, X 3 X5 8 ), (x 4 , £ 24 , X 2 ) A (x 2 X e7 , X 2 V, X 6 X 24 ), 

(x 8 ,x 58 ,a; 5 ,w,A) A ( X5XW}U , XlU}UXlX ^ u V**ll^ jU ;WZl^) } (5.4) 
(x 7 ,x 67 ,x 6 ,v,£,/i, 1/) A 

(. ,„ „ „ „ „ „ , .1/2 X2^7^ , .1/2 ^3^4^ , .Xi^M^A X 2 UV^ 

{ux 4 x 67 ,x 2 v,x 7 x 24 ,uj X g ,u> "^^, w ^^ , xpiT^ 

Substituting formulae (5.1-5.4) in (2.11) one can show that the whole 
W function will be invariant under r and p transformations from the group 
G. The gauge factors before curly brackets in (2.11) are chosen in a such 
way to cancel multipliers coming after the Fourier transformation of the sum 
in (2.11). The action of any element from G on coordinates can be easily 
calculated by a composition of formulae (5.3,5.4). 

Let us briefly comment on these formulae. It is easy to see that after the 
p transformation, the modulus parameter a (see 3.9) transforms as: a — ► o -1 . 
So it is not invariant under elements from the group G containing an odd 
number of p transformations. In some sense this is a pseudo-invariant (it is 
obvious that any combination like a + a" 1 will be invariant under all transfor- 
mations from G) . Introducing "angles" as arguments of the elliptic functions 
(see Appendix) we can define a "crossing" transformation which in fact cor- 
responds to some spin permutation like {a\efg\bcd\h} — > {a\fge\cdb\h}) and 
etc. 

In the next section we will show that the weight functions of the lattice 
should have alternating values of moduli: a and a -1 in all directions in order 
to provide a commuting family of two-layer transfer-matrices. 



6. Two- layer model 

In this section we will construct a commuting family of two-layer transfer- 
matrices from solutions of the modified tetrahedron equations (2.1). More 
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explicitly, we will show that a composite weight function consisting of eight 
elementary W satisfying (2.1) can be chosen in a such way that it will satisfy 
the tetrahedron equations of mixed type with spin variables placed in the 
corner sites, in the middles of edges and in the centers of all faces of the 
composite cube. 

Consider a composite weight W which consists of eight independent weight 
functions W a , W e , Wf, W g , W b , W c , W d , W h (see Fig. 6.1) and the same for 
W, W and W". 




Fig. 6.1 A composite weight W. 
Let us demand that four weights W, W, W" and W" should satisfy to 
usual tetrahedron equations of mixed type: 



W W W W" = W" W W W. 



(6.1) 



So constructing from W and W" two-layer transfer matrices T(W) and T(W') 
in a usual way we come to a commuting relation: 



[T(W),T(W)\ = 



(6.2) 



Then one can show that it is equivalent to the validity of the following 
16 equations of type (2.1): 



W h W' c W!W™ = W h W c W e W a , W a W f W b W h = WaW'fWZWil', 

WcW'cW^w'j' = w c w'^" c w"g, WfWfWfWZ = WfWfW/ws', 
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w b w>w>W" = w h w' g W" a w^ w e w' d w" h wn' = w e w' d w£W'i 

WdW'Jv'lW,! 1 = W d W' e Wi'W", W g W' h WHW'; = WgWlW'lWl", (6.3) 
W h W' h W^W'f = W h W h W d W" f \ W a W' a W" g W" = W a W' a W' g 'W'^ 
W e W' e W'cW{!' = W e WlW?W", WbWlW'/W'f = W b w' h w" f W'f, 

WfW'WlW" = W f W' a w:w>", W c W' h w^wi" = w c w h w" h w'l 

W g W' g W^W' g = W g W' g W" g W' g \ W d W d W d W%> = W d W' d W' d W d . 

An arrangement of spin variables in relations (6.1) is the same as in the 
modified tetrahedron equations (2.1). Note that we must introduce a set 
of additional weights W a , . . . , Wh which appear only in intermediate steps 
when proving relations (6.1) and cancel after 16-fold application of relations 
(6.3). 

Not let us analyze relations (6.3). The simplest way to satisfy them is to 
put W g = W e = W f = W h = W &nd_W d = W b = W c = W a =W. In this 
case we must demand that W and W should satisfy relations (2.1) and its 
dual version. This is the case considered in paper [14] where it was shown 
that then we have only two solutions: namely the Bazhanov-Baxter model 
and the "static limit" elliptic model. So we should try to avoid this situation. 

But using results of the Section 3 we can directly analyze solutions of 
system (6.3). To simplify our discussion we will omit bellow all gauge factors 
before curly brackets entering in formulae for weight function (2.11) and re- 
lations for corresponding variables u, v, A, £, /i, v. They can be easily restored 
by using relations like (3.19-3.21). 

In Section 3 we have seen that all relations containing the coordinates 
x follow (up to some roots of unity) from relations (3.11,3.22). So we will 
analyze only those constraints which contain variables x ai , x ai j, . . . , x^, x^j. 
Introducing notations like (3.9-3.10) for the rapidity variables for the letters: 
a, ... ,h we come to the following relations: 

Jcl = 1g\i 3h\ = Hi, 3e\ = ^a\i 3d\ = % } \i 3b2 = ^92' -?/2 = ^ a2 > 
jh2 = "ic2i jd2 = "ie2i ja.3 = ig 3 , jf 3 = Hz, je3 = ic3i 3d3 = ^3(6.4) 

a g = a e = df = ah = a d x = a^ 1 = a^ 1 = a^ 1 = a = a' = a" = a'" (6.5) 

and the same for all i',i",i'",j',j",j m . 

Relations (6.5) show that modulus parameters a« associated with elemen- 
tary cubes of the lattice have a checkerboard structure in all directions. 
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We also have 24 relations like (3.11): 



4 = 4 »i 4 = 4 > 4 = 4 " 4 l =4 l" 4 l =4 " 4 n =4 I" 

°91 °92 ' °92 t 92' °93 92' 91 93 ' 93 93' 91 °91 ' 



]d\ — 3d2 i Jd2 ~ Jd2i 3d% — 3d2i 3d\ — 3d3 ; 3d3 — Jd^j 3d\ — 3d\ t 



(6.6) 
in 

(6.7) 



and 

• / • • // • / • II • / • // • / • /// • // • /// In c 

<-c2 — ''c2i *c3 — *c2) *c3 — *c3; ^t>3 — *b3? ^bl — *&3 ? 1-bl — i-bl i \y-t 



• / • // • /// • / • /// 



(6.9) 



(6.10) 



At last, we have 16 tetrahedron constraints like (3.22): 

3^58 x e5 x a2A _ , , x a2A X f6 X b67 — 1 X f% X f %7 X f6 X d\3 
X h Tl 777" — tU; ~~x 7 77T — - 1 -; ~x~„ 7 77 777 

13 ^ess X a2 a2 3-/67 Xfl 13 ' 67 X f6 X Hi Xd l 

II III I III I II III 

, , Xb67 x a24 X a2 _ 1 x eb X d\3 x h\3 _ i 3^ 13 X e5 X e58 X e5 

^ 7 77 77T — J-j x 7 77" — J-; "jT 7 77 777" 

"^924 Xa 2 X a24 e ^ x d\ x h\3 dl 3?e58 -^e5 x e§g 

x hJ3 X d\3 X _J6 — I x e2 x eb8 X c5S _ ]_ ^667 £/67 £/2 _ 

x h\3 x d\ x f & 7 e58 3: e2 X c5g Xb 67 Xf 2 Xf 67 

, , X f&Z X o-2 X a24 X e2 _ 1 x a2 X a24 X c5S _ 1 3^13 Xrfj 3^q!i3 -^l _ 1 

<jJ ~X~f 7 77 777" — J-j x 7 77~~ — J-; "tv^ 7 77 777" — - 1 -) 

72 3? a2 4 3? a2 3? e tjg a24 3J a2 3^ g2 4 rf l X d\3 x d\ x d\3 

II I III I III „ „ // 

^^# = 1, wf^^fF = 1, a; ^13 ^67 = 1 924 924 =L 

X c58 3^g 24 924 3^ 67 c58 X^ 13 :T S24 aj ff24 

Since relations (6.4-6.5) leave 13 degrees of freedom in each of the composite 
weights (12 "spectral" and one modulus parameters), we have 40 relations 
(2.3-2.7) on 12 x 4 + 1 = 49 variables. 

One can show that all other relations constrained x and x are satisfied 
automatically up to some roots of unity provided that relations (6.4-6.10) 
are valid. The proof of this fact is quite tedious and we will not give it 
here. But we must be sure that there exists a consistent choice of degrees 
of u in all parameters so that all 16 equations (6.3) will be valid. Below we 
will give a simplified solution of system (6.4-6.10) and the whole system of 
equations (6.3). It shows we should find a general solution only for relations 
(6.4-6.10) and then all degrees of u in parameters x and x can be restored 
by a continuous deformation of the parameters for the simplified solution. 
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Now let us turn to system (6.4-6.10). Unfortunately, up to now we did 
not succeed in obtaining a complete understanding of the underlying geomet- 
rical picture (if any). If we put a = 1 then we come to the inhomogeneous 
variant of the TV-state generalization of the Zamolodchikov model proposed 
by Bazhanov and Baxter. In this case we have a natural geometrical picture. 
Each elementary weight is described by the three angles between the three 
planes. Then a two-layer composite cube is defined by six intersecting planes 
and their mutual orientation depends on nine angle variables. All system 
(6.4-6.10) has 13 degrees of freedom and therefore we have a four-parametric 
commuting family of two-layer transfer-matrices. 

Unfortunately, the case a^ 1 looks much more complicated due to the 
absence of a geometrical analogy. Nevertheless, we have analyzed it quite 
carefully. The results can be formulated as follows: system (6.4-6.10) has 
12 degrees of freedom (11 "angle" parameters and one modulus). A com- 
posite two-layer cube depends from only from eight independent "angles" 
and modulus. So we have an additional constraint among nine "angles" and 
this constraint will disappear in the limit a — ► 1. We also obtain only a 
three-parametric family of commuting two-layer transfer-matrices. We will 
not write here these formulae since they are rather cumbersome and therefore 
practically useless. Nevertheless, we are sure that there should exist a ho- 
mogeneous solution to the system (6. 4-6. 10), so that all weights W a , . . . , Wh 
depend on the same three spectral parameters and modulus) which will gen- 
eralize the model of Bazhanov-Baxter and "static limit" elliptic model of 
[14]. 

Now we will consider another limit case of system (6.4-6.10). Demand 
that all parameters x a , . . . , Xh are expressed in terms of x g by the following 
formulae: 

x al = ^l x g\i x a2 = UJ / X g2> X a5 = ^ l X g^i X aQ = *-/ X gQi (6.11) 
„. _ , ,1/2 X 913 ~ - , ,-1/2 ^£24 ry. — , ,1/2 X 958 ~ _ , ,-1/2 ^967 



^91 ' x 92 x 95 x 96 



X 



el 



1/XgQ, X e2 — UJ/Xg 5 , X e5 — UJ/Xg 2 , Xf, g — 1/Xg-y, (6 .12) 



„ — , ,l/2 X g67 „ — ,,,1/2 X 958 „ _, ,1/2^924 „ _ , ,-1/2 X 9 13 

x el3 ~~ ^ ~r ' X e24 — <" ~r ' X e58 ~~ ^ "r 1 x e67 ~ <*> ~^ j 

x 96 x 95 x 92 x 91 



X fl = Uj/Xg 5 , X f2 = 1/Xg 6 , Xf 5 = U/Xg v X/g = l/Xg^ (6.13) 

r _ 1/2 ^958 r . - , , 1/2 X 967 T . _ 1/2 ^913 r . - UJ -l/2 X 924 

X /l3 M X g5 ' X /24 M Xg Q ' X /58 M X gi ' X /67 ^ X g2 ' 
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Xd\ — 1/^92' Xd 2 ~ *-/ X 9li X d5 — *-/ x gQi x dQ — ^-/ x g^i (6.14) 

rr , — , ,-1/2 ^£24 ry. , ._ , ,1/2 X <?13 ~ , _ , ,l/2 X ff67 ™, _ , ,-l/2 X ff58 
^dl3 — u ~~r 1 x d24 — W ~r 1 x d58 — w ~f 5 x d67 ~ W ~r 5 

Jjg 2 U,g t J^g 6 Xg 5 



Xb\ — Xg 5 /U, Xb2 — x g$, x bb ~ x g\i x b% — x g2l UJ i (6.15) 

^613 = x g^l UJ i x b2A = UJX g§Ti x b$& — x g\^i x b67 = X 924/ UJ ' 

x c\ x gQi x c2 x g§l^i x c5 x g2i x c6 x g\l^i ^D.IDJ 

el3 -^g^Ti • v c2A 958' a 'c58 * / '<?24' u, c67 • h g\y,l i 

x h\ = x g2l UJ i X h2 = x g\l UJ i x h5 = x g§i x h% = x gz > l UJ i (6.17) 

x hl3 = X g2Al UJ ' X h2A = x gi%j x h58 = X g§-ji X h67 = X g5S/ UJ 

and same for x', x", x"' and x x', x", x'" . Then all relations (6.4-6.10) will be 
satisfied automatically as a consequence of the relations for x g (6.6) and the 
tetrahedron constraints for x g (see (3.22)). Moreover, all relations containing 
x will be satisfied also. So we can use formulae (2.8-3.4) in order to determine 
all u multipliers. Now we will show that transformations (6.11-3.4) have a 
natural interpretation. 

First consider a transformation of x given by (6.14). Let us denote a 
weight function from this set of parameters as W^\ Then one can show that 
W and W^ 1 ' satisfy to the following inversion relation: 

£ W{a\efg\bcd\a)W {i \a\bcd\efg\h) = ®(x)8 a , h , (6.18) 

where <&(x) is some scalar factor which can be calculated. 
Define also the following automorphism 9 on x's as: 

0(^1, ^2,^5, x 6, X 13, x 2Ai x h8i ^67 ) = 

(xi/u, x 2 /uj, x 5 /u, x 6 /u, xi 3 /u 2 , x 24 , x 5 s/v, x 67 /u). (6.19) 

Then all relations (6.11-6.17) can be rewritten in the following compact form: 

W a = wt\ W h = W®, W c = wf\ W d = W®, (6.20) 
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where the action of the transformations r and p on x's was defined in Section 
5 and the action of their compositions can be easily calculated. We call the 
choice of weights in the form (6.20) as the "inverse" variant of two-layer 
model. 

Note that we were forced to introduce automorphism 6 in order to satisfy 
all the systems of relations for x's and x's coming from (6.3). This is a very 
overdetermined system of scalar equations and it is remarkable that we have 
any solution at all. 

7. Conclusion 

In this paper we have found a new solution to the modified tetrahedron 
equations in the framework of the so called "Body-Centered-Cube" ansatz. 
The weight functions depend on three angle-like parameter $i and one elliptic 
modulus k. We have also constructed the elliptic parameterization of the 
weight function. Note that in the limit k ^ Owe obtain the Baxter-Bazhanov 
model . If the "static limit" condition 61 + 62 + 63 = 2K with K being 
the full elliptic integral is satisfied than we obtain the model proposed in 
[9,14]. As was mentioned above we were forced to achieve the integrability 
of the model without the validity of "dual" version of MTE to be able to 
generalize the Baxter-Bazhanov model and not to restrict ourself onto "static- 
limit" condition. The constructed integrable two-layer model has the spin 
variables in the cube vertices, the middle points of the cube edges and the 
centers of the cube faces as was shown in fig. 6.1 . So, this model is of 
mixed type. The weight function consist of eight weights W a , Wb, ...Wh which 
must satisfy the system of the sixteen sets of equations (6.3). The solution 
to these equations appeared to have twelve free parameters (including the 
modulus k) but the formulae are rather cumbersome and we do not write 
them. We have also considered one particular case when all four pairs of 
the opposite weights (a-h,b-e,c-f,d-g) are connected with each other by the 
inversion. Unfortunately, we did not succeed to find a geometrical picture like 
in the case of the Baxter-Bazhanov model. Nevertheless some of the formulae 
presented in Appendix resemble the formulae for the spherical triangle taking 
place for the angles 61,82,63 in the Baxter-Bazhanov model. So, we hope 
that some appropriate geometrical picture does exist. We also hope that 
besides the "inverse" model some another two-layer "homogeneous" model 
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which might have a rich physical content can be found. We expect to do this 

elsewhere. 
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8. Appendix 

The formulae for "modified spherical triangle" are more aesthetic if we intro- 
duce a parameterization for Ti,T 2 ,T 3 in terms of the elliptic functions of the 
modulus k — (1 — m)/(l + m). Besides, some of the formulae bellow resem- 
ble corresponding formulae in the spherical trigonometry. We shall use the 
ordinary Jacobi's elliptic functions sn, en and dn, theta - and eta - functions 
H , . Let us also introduce the following notations: 

/ , x , / sn(u, k) .„ „ . 

tn(u,k) = k' \' (8.1) 

cn(w, k)dn{u, k) 

and 

2k 1/2 — 2ik 1/2 

tn L (w, k) = tn((l + k)u, — — -), tn L (w, k) = tn((l - k)u, — ) 

tn L (u, k)tn L (u, k) = tn 2 (w, k). (8.2) 

The last two functions represent the Landen's transformation r — ► r/2 and 
t — > r/2 + 1. The parameters k and m we regard to be fixed forever and 
so we shall omit them everywhere. Also we shall use the indexes p, q, r for 
1,2,3 

Defining 9 r so as 

T r = tn(0 r /2), (8.3) 

then we obtain 

S r = sn(9 r ), C r = cn(0 r ), D r = dn(9 r ). (8.4) 

Having got the angles r , define the spherical excesses as usual: 

6»i + 6 2 + 3 
«o = = A-, « r = u r — ckq, (8.5) 
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where /C is the complete elliptic integral of the first kind of the modulus k. 
It can be shown that the following expressions are valid: 



, ^ ,., H(a )H(a 1 )H(a 2 )H(a 3 ) 



sin ao/2 



9(0)9(fl 1 )9(fl 2 )9(# 3 
cos 2 a /2 



,, M 9(a )9(a 1 )9(« 2 )9(a 3 ) 



9(0)9(^)9(^)9(^3) ' 
tan— = A;ysn(o;o)sn(ai)sn(Q;2)sn(Q;3); (8.6) 



., , H(a )H(a r )Q(a p )&(a q ) 



sin a r /2 



Q(0)Q(6 r )H(6 p )H(6 q 
cos 2 a r /2 



,, Q(a )Q(a r )H(a p )H(a q ) 



e(0)e(9 r )H(9 p )H(9 q 



a r 
tan- = . 



sn(a )sn(Q; r 



sn^a p jsn^a g j 



17) 



Similarly one can obtain for the linear excesses: 



2 Pq _ tn L (a 1 /2)tn L (a 2 /2)tn L (a 3 /2) 



tn L (a /2) 
tan 2 — = tn L( a o/ 2 ) tn L(a P /'2)tn L (a q /2) 

2 ~ tn L (a r /2) ' [ " J 

The formulae for the /3 can be obtained from (8.8) by replacing all the tn^ by 
tnj, (it is the common rule: replacing k by —k is equivalent to the modular 
transformation r to r + 2). We complete our list of the formulae by 

^§ = tn L (V2)tn L (^/2), 
sin fjQ 

sm p sin p r "" sin p sin p r 

Now let us give the explicit form of the Boltzmann weight function for the 
case N = 2. The calculations resembles ones from the Appendix of [14] and 
we give here only the answer. In our previous sections all the spins belong 
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to Zm and take the values 0, 1 N — 1. For the case N = 2 it is convenient 



/, J-, .. ., 



to deal with multiplicative spins (— l) a = ±1 instead of a = 0, 1. Further we 
will imply the multiplicative spins as the arguments of the Boltzmann weight 
functions. 

Introduce for shortness some conventional notations: 



and 



y/snK/2), a = yJcd(ai/2), ? = 0,1,2,3; (8.10) 



Pq — C0C1C2C3, P r — CQC r s p s q ; 

Qo — S0S1S2SS, Q r = SoS r CpC q ; 

Ri = SiCi, Si = k , 1 = 0,1,2,3. (8.11) 

Then for the weight W given by (2.11) there exists the following identity: 

1 



W 



-W : 



Pq — Qo 
and Wspqr is given by the following table: 



SPQR 



8.12) 



abeh 


acfh 


adgh 


W S PQ R (a\efg\bcd\h) 


+ 


+ 


+ 


P — abcdQ 


— 


+ 


+ 


Ri + abcdSi 


+ 


— 


+ 


R 2 + abcdS 2 


+ 


+ 


— 


i?3 + abcdSs 


+ 


— 


— 


abPi + cdQi 


— 


+ 


— 


acP 2 + bdQ 2 


— 


— 


+ 


adP 3 + bcQ 3 


— 


— 


— 


R + abcdS 



Note that for Wspqr all the symmetry properties are quite evident. 
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